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The effect of low energy excitation on the dynamics of Josephson phase in SNS junc- 
tion is investigated. From the microscopic Hamiltonian, the effective action for the 
phase variable is derived. The retardation effects due to low energy excitation in nor- 
mal(N) region are seen to play important roles in the dynamics of the phase. By the self 
consistent harmonic approximation, renormalization of mass and dissipation constant 
are calculated, revealing the enhancement of the former and the suppression of the 
latter in general. Various situation appears depending on these renormalized values. 
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§1. Introduction 

The quantum mechanics of the phase of a superconductor in Josephson junction has been 
studied both theoretically and experimentally by many authors.^ This problem was first 
treated phenomenologically in the pioneering work of Caldeira and Leggett based on the 
functional integral approach.!' i-* Since the functional integral approach gives a clear view and 
the insight to the physics, it was developed thereafter to give the microscopic foundation to 
the quantum mechanics of the Josephson phase.BBB In these studies, the dynamics of the 
phase variable is described by the effective action, which can be derived from microscopic 
Hamiltonian by eliminating the electron field. It was shown that the effective action consists 
of the mass term which describes the charging energy, E c , of the junction and the terms 
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which express the effect of tunneling of quasi-particle and Cooper pair across the junction. 
The latter two terms, which express the effect of tunneling, can be viewed as the interaction 
between the phase at different imaginary times. 

In the case of superconductor-insulator-superconductor(SIS) junction with large capaci- 
tance, the mass of the phase is large and thus the time scale for the motion of the phase is 
long compared to that of the interaction term. Therefore the "adiabatic" approximation can 
be justified and the kernel expressing the interaction between the phase at different imagi- 
nary time can be approximated by 5-function. With this approximation, the action reduces 
to that of a particle in a cosine potential with Ohmic dissipation. In this case, macroscopic 
quantum tunneling of the phase is predicted and experimentally observed.ll 3 i'iS 

For a SIS junction with small capacitance, Geigenmiiller and Ueda0) investigated the 
non-adiabatic effect on the effective critical current of the junction by use of self consistent 
harmonic approximation(SCHA). They compared the results of SCHA with those of adia- 
batic approximation and find the qualitative agreement between the results obtained by the 
two approximations even if the charging energy is larger than the gap of the superconduc- 
tor. Quantitatively, however, the difference between the adiabatic approximation and SCHA 
turns out to be quite large. 

As for the superconductor-normal metal-superconductor (SNS) junction, it is now possible 
to fabricate a junction with a capacitance small enough to observe macroscopic quantum 
effect.© In contrast to SIS junction, low energy excitation in normal(N) region gives rise to 
the long time behavior of the kernel, K(t), and hence K(t) may not always be approximated 
by a ^-function. 

The classical limit of the SNS junction with two dimensional electron gas(2DEG) with 
diffusive motion was investigated by Kresin,© who calculated the critical current, j c . In his 
approach, the charging energy, E c , of the junction is assumed to be zero and, correspondingly, 
the phase valuable is independent of the imaginary time, i.e. 9{r) = Oq. In this case, 
Josephson coupling energy is given by Ej = §{f &tK(t), and the critical current, j c , is 
obtained as j c = 2eEj. When the charging energy becomes finite, however, the quantum 
fluctuation of phase must be taken into account and we expect the retardation effect of K{r) 
to appear in the dynamics of the phase. In this paper, we will investigate this retardation 
effect on the tunneling rate of the phase in the neighborhood of critical current, where the 
tunneling barrier height becomes small. 

Apart from the problem of Josephson junction, the tunneling rate in the small barrier 
height region was studied by Kleinert when there is no dissipation and the retardation 
effect.0 In the present paper, we will extend his method to the case with the dissipation 
and retardation effect to calculate the tunneling rate in the neighborhood of the critical 



2 



current. 

To investigate how the retardation of the kernel affects the decay rate of the metastable 
state, we first derive the effective action for the Josephson phase in SNS junction from the 
microscopic Hamiltonian, following the procedure outlined by Ambegaokar et a/.&&. Then 
we apply self consistent harmonic approximation to the action, and calculate the renormal- 
ization of the mass, dissipation constant and attempt frequency. With these renormalized 
parameters, we will study the tunneling rate in the neighborhood of the critical current. 

The organization of the paper is as follows. In the next section, we will derive the effective 
action for Josephson phase in SNS junction in the perturbation theory with respect to the 
mixing matrix element between S and N region. In section 3, we will explain the self- 
consistent harmonic approximation applied to the effective action derived in section 2. The 
results and discussion will be given in section 4 and a brief summary in section 5. 

We take units, h = k B = 1, unless noted. 

§2. Effective Action for SNS Junction 

In this section, we will derive the effective action from microscopic Hamiltonian following 
Ambegaokar et al® 

2. 1 Derivation of the Effective Action 

The model of the junction is as follows(See Fig. 1): Two superconducting regions(S region), 
which will be referred as L- and .R-region, are connected by a normal region(N region) of 
length, L, and width, W. The S regions are assumed to be clean bulk superconductor and N 
region to be 2DEG with impurities. Finite transfer integral, t rr i, of the electrons is assumed 
through the boundaries of S and N regions. The transfer integral can be taken as real without 
a loss of generality. We will take z axes perpendicular to the 2DEG. Other axes, rn and r± 
are taken to be parallel to and perpendicular to the junction direction, respectively, as shown 
in Fig. 1. 

The Hamiltonian of the system is given as follows. 
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Fig. 1. A schematic picture of SNS junction. The top view(a) and the cross sectional view(b) are shown. 
The superconducting(S) regions, L and R, are connected to 2DEG of normal(N) region by the transfer 
integral t rr >. The width and the length of the junction are W and L, respectively. The 2DEG plane lies 
at z = 0. 



where g, m e , [i and — e are interaction constant, effective mass of an electron, chemical 
potential and electron charge, respectively. The electron field and scalar potential are ex- 
pressed by ip(r) and <f>(r), respectively. The impurities in N region is expressed by a potential 



Here JVj, V- nnp and rj are the number of impurities, the strength of the impurity potential 
and the position of the j-th impurity, respectively. Since the wave function of the electrons 
in N-region is restricted to a narrow region around z = as in Fig. lb, the scalar potential 
in N-region is evaluated at z — 0. In eq.(2.1), the effect of vector potential is not explicitly 
taken into account, because it is not necessary in the following discussions. 
The action of the system is defined by, 
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In eq.(2.2), the Wick rotation is performed, i.e. it — > r and — > i<p. 
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First, we eliminate the quartic interaction term in Hamiltonian, eq.(2.1), by the Hubbard- 
Stratonovich transformation, 



1 = J VA*VA exp -±- J dx (A(x) - g^ L (x)^(x)) (A*(x) - g^W^] 
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where x denotes both imaginary time, r, and position, r. With this transformation, the 
action becomes, 
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The Green's function, G ( 6x6 matrix spanned both in the space of L, N, R and 

Nambu space. The electron fields are expressed as, 

(x) 

* (x) = 
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The inverse of the Green's function for electron system, G~ l (x,x'), is given by. 
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where n, r 3 are the Pauli matrices and 1 is a unit matrix and #L(i?) is the phase of the 
superconducting order parameter, ^l(r) = Ac 



L(R) 



Before tracing out the electron degrees of freedom from the action, eq.(2.4), we perform a 
transformation, 

*L{N,R) (X) = S^'V^VlWR) {x) , (2.8) 

to express the phase explicitly and then eliminate them from the order parameters in eq.(2.7). 
Here we have also introduced a phase in normal region which satisfies©, 

d 



Or 



9i<i(x, r) = 2e<j) (x, r) 



(2.9) 



5 



to get rid of the scalar potential in N- region. Note that the scalar potential in eq^^) is 
considered as <f>(x) = <f){r\\,r±,z = 0, r), since the electron gas in N-region is confined in 
the narrow region around z — 0. With eq.( |2.9| ) the phase can be defined continuously from 
S-region to N-region at the boundary. (See eq. ( |2.17| .)) Then we arrive at the action, 



S e i = / dxdx'V* (x) —G (x, x ) V (x' 



(2.10) 



where the inverse of Green's function is expressed as, 
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The diagonal elements are, 
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Tracing out the electron fields from eq.( |2.10|) , we obtain the action as follows to the second 
order of the matrix element of the transfer integral, 
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In eq.(2.13), tr means taking trace with respect to the matrix element of the Green's function 
and integrating with respect to x,x f . 

Now the action, S, depends on A, <p and 6, which are independent fluctuating variables. 
Fortunately, these quantities, except the phase 8, have their mean field values with only 
small fluctuations around them, which we ignore. 

First we investigate the zero-th order term of eq.(2.13). The action for L, N, R region are 
defined by, 
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Expanding the Green's function in L- and R-region with respect to V0(x) and d T 6(x)—ie<p(x) 
to the second order, we obtain as the relevant part of the effective action for S-region, 
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where ^l(r) is the volume of the superconductor, p s is the superffuid density and N is the 
density of states per spin at Fermi energy.^ 

From eq.(2.16), we see that the action is extremal when the phase variable satisfies the 
Josephson relation,© 

^-9 (x) = 2ecj) (x) . (2.17) 
Or 

With eqs.(|2.9|) and (|2.17|) , the phase of the S region and that of the N region are connected 
smoothly at the boundary. 
On the other hand in N-region, by expanding the action with V0(x), we obtain, 
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where Trj^g (x, x') = —(T T j a (r, r)j /3 (r', r')) is current correlation function in N region, p is 
the density of the electrons and indices a, (3 stand for r± and r\\. Here T T is the time ordering 
operator and j a (r,r) is the current operator. Since we have taken the phase which satisfies 
eq. (|2.9| ), the contribution from density correlation function does not appear in eq. fl2.18p . 
Noting that the characteristic decay length of Trji^g (x, x') is the electron mean free path, /, 
and the variation of the scalar potential is small, 9{x) and 9(x') in the first term of eq. ( |2. 18| ) 
can be evaluated at the same point, i.e., r = r'. By performing integration with respect to 
r' in 7T]^ a a(x, x') 7 the effective action in eq.( |2.18|) can be approximated as, 
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where the correlation function, tt(t,t'), is given by, 
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where r tr is the transport scattering time of the electron due to elastic scattering and uT tr <C 1 
is assumed. Weak localization effect is not taken into account here. 
In terms of the phase difference, 
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which is assumed to be uniform in r± direction, eq.(2.19) becomes, 
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where a = Rq/Rn- Here Rq = h/(4e ) ~ 6.45kf2 is the quantum resistance and R 
L/{Wa) is the resistance of the junction. In eq.( [2.22|) , we assumed V# (r, r) ~ Q- (r) /L, 
which means the electric field is constant in N region. The term eq.( [2.22| ) has been derived 
by Schon and Zaikin for SNS junction.0 

Next we investigate the effect of the matrix element of the transfer integral in eq.(2.13). 
The diagrams expressing these processes are listed in Fig. 2 



Fig. 2. The diagrams which appears in the expansion of eq.(2.13) with respect to the matrix element of 
the transfer integral t. 



The lowest order contribution (See Fig. 2a), is given as, 
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Here we made use of the fact that the matrix element of the transfer integral, t xx >, differs 
from only for x ~ x' and in the neighborhood of the S-N boundary, i.e. t xx > = t8(r\\ — 
r '||)^( r ± — r '±)^( r || ~ r \\)&{ z )- Here r|| is the position of boundary between R(L) and N-region. 
Therefore the integral in eq.(2.23) is taken over r and rj_. The Green's function g^ and g^ 
in eq.(2.23) are the 1,1- and 2,2- component of Nambu Green's function, respectively. Using 



eqs.(2.7a), (gTg) and (pT7|) , we find, 

^L(R) ( X > X ) = — 



5 (x — x') . 



(2.24) 



Due to the existence of the gap, A, in the excitation spectrum, the Green's function becomes 
short ranged in time and space. Thus we approximate it as, G(x — x') oc 5 (x — x'). 

We also use the quasi-classical approximation to the Green's function in N region, which 
leads to, 
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exp 
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\2dr 
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where da? • (a, v) = adr + v ■ dr and the path C is shown in Fig. 3a. Since Gl (x, x') oc 
6 (x — x'), the phase contribution in eq. (|2.25| ) disappears. Thus the term shown in Fig. 2a 
is independent of the phase. 



Fig. 3. Contours for the integrals in (a) eq.(2.25), (b) eq.(2.27) and (c) eq.(2.30b). 



Next we turn to the fourth order diagrams, i.e. Fig. 2b and Fig. 2c. The first diagram is 
expressed as, 

t 4 r 

Sa = — J dx 1 dx 2 dx 3 dx 4 gL]- (zi, x 2 ) g N ^ (x 2 , x 3 ) g L] - (x 3 , x 4 ) g m (x 4 , xi) 
+hi (xi,x 2 ) g Nl (x 2 , x 3 ) g Ll (x 3 , x 4 ) g N[ (x 4 , x x ) 
+h (xi,x 2 ) g Nl (x 2 , x 3 ) fl (x 3 , x 4 ) g N ] {x 4 , x x ) 

+fl (xi, x 2 ) g Ni (x 2 , x 3 ) f L (x 3 , x 4 ) g m (x 4 , x x ) . (2.26) 

Here /, /* denote 1,2- and 2,1- component of the Green's function. The first and the sec- 
ond term do not contain the phase variable, because the contributions from g^(x,x') and 
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g^(x',x) cancels each other. In the third and the forth term, the phase factor is given by, 

' 1 89 1 . 
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with the path shown in Fig. 3b. In principle, the contour, C, for the integrand must be 
determined from the extremal path for the electron. Here we approximate it with the line 
connecting the terminal points, x, x' as shown in Fig.3b. Noting the Josephson relation, 
eq.( |2.9| ), and that the junction is uniform in r± direction, the phase contribution in eq. ( [2.271 ) 
will vanish. This reflects the assumption we have made that there is no voltage drop across 
the S-N boundary. So the diagram in Fig. 2b can also be neglected. 

The last diagram, Fig. 2c, is the only relevant term in our analysis. From this, we obtain 
the action as 

t 4 
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+f L (x!,x 2 ) g N i (x 2 , x 3 ) f* R (x 3 , x^ g m (x 4 , an) 
+ft (xi, x 2 ) g Nl (x 2 , x 3 ) f R (x 3 , x 4 ) g m (x 4 , x x ) 
+ [L^R}. 



(2.28) 



In eq.(2.28), we find the third and the forth terms, which express Cooper pair transfer 
through the N region, are relevant and lead to, 



Sj = t A J dx!dx2dx3dx 4 |/£ (x x , x 2 ) g Ni (x 2 , x 3 ) f R (x 3 , x 4 ) g m (x 4 , xi) 
+ f* R {xi, x 2 ) g Nl (x 2 , x 3 ) f L (x 3 , x 4 ) g m (x 4 , x x ) . 
For low energy, we obtain, 

Sj = - J dr ± dr' ± dTdT'K(x, x') cos [W c (x, x')} , 
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(2.30a) 
(2.30b) 
(2.30c) 



where the contour, C, is taken as in Fig. 3c and (• • -)i mp denotes impurity average. The 
integral region in eqs. (2.30a) and (2.30c) are on the SN boundary. Now we consider the 
case when the time dependence of the electric field is weak. In this case, the contour can be 
approximated by the average of C3 and C4 of Fig. 3c, since the contribution from path CI 
and C2 in Fig. 3c cancels each other. Thus we obtain, 

6-(t)+6-(t') 



W. 
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(2.31) 
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With this approximation, the action becomes, 



rP 

Sj = — drdr K(r — r) cos 
Jo 
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(2.32) 



where Kir — r') = J dr±dr' ± K(x,x'). From eqs. Q2.22] ) and ( |2.32| ), we obtain contributions 
to the action from the electron field as, 



Sei = ^E dn^9Z n - f drdr'K (r - r') cos 
Zp „ Z7T Jo 



9- (r) + 9- (r f ) 
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With the help of eq.( |2.9| ), the contribution from S em , eq.(2.4c), is given as, 



3 C (d9-\ 2 

where C is the capacitance of the junction. In the derivation of eq. ( 2.34ft , we have assumed 
that 9 is uniform in r± direction and used the approximation, V# ~ 6~ /L. The total action 
for the phase is now reduced to the effective action, S e s, which is written with the single 
degrees of freedom, 9, as, 

'9(t)+9(t')' 



^ = ^"E^g^-n- 1 drdr'K ■ 



T COS 



#(7(2,35a) 



90,1 = ra^n + T^KI, (2.35b) 
m = £ = (2.35c) 
a = (2.35d) 

it AT 

where we have denoted 9~ as 9. The last term of eq. (2.35a) is added to the action to express 
the bias current, j. 

When the kernel K(r) is a delta function, the effective action, eq.(2.35) reduces to that of 
the phenomenological model of resistively and capacitively shunted junction(RCSJ), which 
is a good description for SNS junction in the classical limit. 

The action, eq. (2.35a), differs from that of Josephson junction in two respect. First, there is 
an Ohmic dissipation term, a, in eq. (2.35b), whereas in the action of Josephson junction the 
dissipation term is given by sinusoidal form reflecting the discrete transfer of charge between 
the electrodes. Second, the kernels are short ranged in Josephson junction, whereas, in SNS 
junction, the kernel, K(r), has a long time tail due to the low energy excitation of the N 
region. 

The effective action, eq.(2.35), is derived for the first time in SNS junction. We will 
compare this action to the similar system of SINIS junction, which is derived by Guinea and 
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Schon.El'Ill* Due to the fact that there is a voltage drop across the SIN boundaries, there are 
two degrees of freedom, 9r and 9l, which represents the phase difference across right and 
left SIN boundaries, respectively. The contribution of the diagram in Fig. 2b leads to the 
dissipative term in this case. Josephson coupling term, however, leads to 

- / dTdr'Kir - t') cos (6r (t) - 6 L (r')) . (2.36) 
Jo 

Here the kernel Kir — r') is given by the pair propagator as in our results, eq. (2.35a), but 
the difference of the argument of the cosine term in eq.( |2.36j ) is to be noted. 



Finally we comment on the assumption that V# can be approximated as a constant in 
N-region. This means constant electric field is assumed in N-region, although the electric 
field must be determined self-consistently at each point in N-region in principle. Since our 
aim is to investigate the retardation effect of K(t), the essential physics will be seen in this 
approximation. 

In the next section, we will discuss the imaginary part of the free energy based on the 
partition function, 

Z = J V6exp(-S cS ) . (2.37) 

Before going into the discussion of tunneling rate, we will discuss the time dependence of 
the kernel K(t) in the following sub-section. 

2.2 Time Dependence of the Kernel 

The long time behavior of the kernel K(r) is determined by the low energy excitation of the 
system. In Josephson junction, K(t) can be approximated by 5-function due to the existence 
of gap, A, in the excitation spectrum. In SNS junction, however, time dependence of K(t) 
is determined from the low energy excitation of the N region. Various situation is expected 
depending on the difference in the electronic properties. In this paper, we will consider the 
case where the motion of the electrons in 2DEG is diffusive. Actually, the experiments have 
been carried out even in the strong localization regime. We expect, however, the essential 
aspects of physics will be clarified by studying such a diffusive case. 

The kernel K{r) is given in the form, 

K(r) = E J0 k(L,T). (2.38) 

Here k(L, r) is the real space representation of the pair propagator, which is normalized 
appropriately. Since we are interested in the dependence of K(t) on diffusion constant, D, 
the other factors, such as density of states, transfer integral and so on are included in Ejq, 
which is given by 

Em = i*^m, (2 . 39) 

7T L z 
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where N(0) is the density of state of 2DEG per spin. 

In the region where the electron motion is diffusive, the pair propagator, k(q, io> n ), is given 
for ujr tT <C 1 and q <C / _1 as, 

k(q, kj n ) = \Y, — 1 2 , (2.40) 

where v n = (2n + l)irT and i/ c is the cutoff frequency. The diffusion constant is given by 
D = VpT tT /2 = ixa/(e 2 m e ) = 2La/m e W with a = Rq/Rn- The function, k(L,r), is given 
by Fourier transform of eq. ( |2.40| ) with a suitable boundary condition. By the boundary 
condition that the current normal to the SN boundary is zerc00' 1 , we obtain 

1 CO 

k(L,iu n ) = - £ (-l) m k(-m,iu n ). (2.41) 



m=— oo 



From the relation, eq.( |2.41| ), K{t) is given as, 



K(t) = EjAY.^H^u) 

P n 

/OO J, . y-^ — LO I T | 
oo - T —^lmk^(L,u J ), (2.42a) 

k {R \L,u) = k(L,uj + iO + ). (2.42b) 
The real time behavior of the kernel, K(it), is derived from K(t) by the analytic continu- 



ation, i.e., t — > it, as, 



K{\t) = K R (t) +iKi(t), (2.43a) 

K R (t) = E JQ f°° —-S2^Lhnk(L,w), (2.43b) 
J-oc 7r 1 — e p 

= -^f-T^ImMl, W ). (2.43c) 



The time dependence of the kernel, K{ir), is determined by £/L, where £ = J D/2nT 
is the coherence length. The coherence length is related to a as ^/I = \J a / 7xm e W LT . 
To make a direct contact with the experiment of Takayanagi et alj^ we take W = 80/im, 
L = 0.4/im, T = 20mK and m e = 0.026m (InAs) and the Nb electrode is used. Here m 
is the free electron mass. With these parameters, £/L, is given as £/L = 1.3 x 10~ 2 y/a. 
Thus we obtain £/L = 0.41 for a = 10 and £/L = 1.3 for a = 100, respectively. The 
time dependence of K{it) for these values are shown in Fig. 4. Here the solid line expresses 
Kr and the dashed line Kj. The line with square represents the kernel and with the circle 
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Fig. 4. The kernel Ku(t) (solid line) and —K](t) (dashed line) as a function of real time t/ (3 for (a) £/L = 
0.41 and (b) £/L = 1.3. The lines with square is for the boundary condition that the current normal to 
the SN boundary is zero and the line with circle represents its L — > oo limit. 

represents the asymptotic form for L — > oo. We take the cutoff frequency v c as v c ~ A. This 
is consistent with the approximation that the Green's function in S-region is a 5 function. 
From the graph, we see little difference for the time dependence between asymptotic form 
and the exact one even for £/L = 1.3. Thus we will use the limiting form hereafter, for the 
computational simplicity. 

§3. Calculation of Tunneling Rate 

In this section, we will explain the method to calculate the tunneling rate of the Josephson 
phase from the metastable state. In our action, due to the time dependence of the Kernel, 
K(r), the potential term and the dissipation term cannot be divided easily. Thus we employ 
the self consistent harmonic approximation as Geigenmiiller and Ueda.0* Applying self- 
consistent harmonic approximation, we first calculate the renormalization of the tunneling 
potential, the mass and the strength of dissipation. With these parameters, the decay rate 
formula for the cosine potential is applied to calculate the tunneling rate. 

In the next sub-section, we will derive the self-consistent equation. 

3. 1 Self Consistent Harmonic Approximation 

In order to separate the potential term and dissipation term in our effective action, we 
use the self consistent harmonic approximation. This will be done by using the variational 
principle, which was developed by Feynman and Kleinert.i^ 
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The variational principle for the free energy is given by. 

F<F tI + ^(S-S tT ) tT , (3.1) 

e _/3Ftr = j V6e- St \ (3.2) 

Here (• • -) tr denotes the average with respect to trial action S tT . This follows from the 
convexity of the Boltzmann weight function with respect to the action. 

To determine the renormalization of potential, mass and strength of dissipation, the direct 
application of eq.(3.2) is not possible. Thus we will define the effective classical potential, 
W(9), and consider the motion of the phase in W(8). The effective classical potential is 
defined asjill 

exp(-(3W(9)) = J VQ{t) exp (S eS ) , (3.3) 

where 0(r) is the fluctuation around the average path, 9 = Jq dr9 (r) which is defined 
by 

6(r) = 9{r) - 9, 

= \ E ^ T e n . (3.4) 

Here Matsubara frequency is given by u n = 2nn/p. With the effective classical potential, 
the partition function is expressed as, 

d0e"W>. (3.5) 

This form resembles that of the classical partition function with potential W. This is the 
reason we call W, the effective classical potential. This is the exact limiting form of the 
partition function for high temperature, as the deviation of the path from 9 is small in this 
case. 

In a damped harmonic oscillator, we can perform the path integral in eq. ( |3.3[ ) exactly, 
since all the integrals are quadratic in the variable 9 n . In our case, however, the path 
integration cannot be carried out, and some kind of approximation is needed. Here we apply 
the variational principle for W. Since exp(— (3W) is convex as a function of W, the same 
discussion as the derivation of variational principle for the free energy can be applied. The 
variational principle for W is given by, 

W < WV + ^(S - StrV = W, (3.6) 

where W tT ' is a trial effective classical potential evaluated with S tr '-Here prime denotes that 
9 degree of freedom is excluded. 
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is real. However, when 



Above discussion holds, as far as the function, exp [-pw (e) 
the bias current j comes closer to the critical current, the variational parameters acquire 
imaginary parts due to the fact that the potential barrier height becomes comparable to 
the zero point energy of the quasi ground state in the potential well. In such situation, the 
inequality, eq.(^.6|), no longer holds. Nevertheless, based on the discussion of Kleinert, we 
expect that the best approximation to the effective classical potential is obtained when the 
derivative of W with respect to (S — S tr ') equals to zero.0* The reason is as follows. A full 
perturbation expansion of exp(— (3W) with respect to SV is given as, 



w{e) = -hnz t 



(3.7) 



/T v " 2(3 

Here (■ • -)tr',c denotes the cumulant average. By definition, above expansion is clearly 
independent of the trial action, S% T >, if we sum the series to infinite order. This independence, 
however, is lost when we truncate the series to the finite order. But we expect the best trial 
action when the approximation shows minimal dependence on S tT >- Thus even in the region 
where the parameters have imaginary parts, we shall use the same variational free energy as 
in the region where the parameters have no imaginary parts. 

Now let us apply eq.( |3.6| ) to the effective action of the SNS junction given by eq.(2.35). 
We take the trial action of the form, 

1 



5tr 



2(3 



E W*- 



n^O 



= m ven uj n 4 



a r 



2tt 



\Ur, 



-\- T7l ren f2 ren 



(3.8a) 



(3.8b) 



With this trial action, we obtain for W, 



H>=-io g n 



i 



2(3 ^ 1 - 



where g(r) 
obtain, 

SW 



—6 - [ dTK(T)e-^ {9{0)+9{T)} cos 9, (3.9) 
2e Jo 



n =i rrw-gn ^ \9on 
z J2n e ~ luJnT 9n/ P- Taking the functional derivative of W with respect to g n , we 



ft E S 9n 
n=l 

0. 



1 1 1 r 13 / 
- + — + -/ drK(r (l 

9n 90n 2 JO V 



(3.10) 



Next, to find out the potential minima,^, we differentiate W with respect to 9, which leads 
to, 



SW ft 3 



56 



[ drK (r) e -ito(o) + *(r)} sin g _ L _ i. £ ^ g ? - g £ _ I j> K (r) c -ifr<o)-*M> cog - Q 
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drK (t) e ± 



I{9(0)+S(r)} 



sin 9 



J_ 

2e 



0. 



(3.11] 



Here we have used the variational condition, eq.(3.10), in the second equality. Above 
conditions, eqs.(3.10) and (3.11), lead to a set of self-consistent equations, 



9 n 9on 

J_ 

2e 



drK(r)(l + e- ibJnT )e 



•|{9(0)+S(r)} 



COS# , 



drK (t) e 



-|{ff(0)+3(r)} 



sin 9 . 



(3.12a) 
(3.12b) 



These equations are formally the same as those derived by Geigenmiiller and Ueda in 
Josephson junction, when we replace K[t) by corresponding terms in the action which has 
been obtained by Ambegaokar et al.0 However, their approach differs from ours in the 
following point: They took 6 n as an independent variable, whereas in our approach j is the 
independent variable. For this reason, they can not discuss the region where the imaginary 
parts appear in the renormalized parameters. Furthermore, by noticing that the bias current 
j is the controllable variable in experiments, our choice seems more natural. 

Our aim is to calculate the tunneling rate of the phase. For that purpose, we have to 
determine the renormalized mass, m ren , renormalized strength of the dissipation, a ren , and 
renormalized attempt frequency, f2 ren , which are defined at low energy. To get these variables, 
we have to perform analytic continuation of the Green's function in eq.(3.12) from imaginary 
frequency to real frequency. We first transform the variables by r = iz and change the 
integration contour as shown in Fig. 5. Next we perform analytic continuation, iuo n — > o;+i0 + , 
and expand the resulting equation with respect to to. Thus we obtain the following four sets 
of equations as shown in the appendix,© 



m r 



m-\ — e 
4 



i POO -i 

3f (0) / dzz 2 e i 9R{z) 



K R (z) sin 



a 



2tt 



1 

2 6 



ni 



ren" r en 



J_ 

2e 



e -|fl(0) 



dze 



dzze 



j9r( z ) 



K 



r{z) sin 



K R (z) sin 
9i (z 



9i [z) 
4 

9i(z) 
4 



— Kj (z) cos 



H dze~* 9R{z) 
Jo 



K R (z) sin 



4 

9i (z) 
4 



— Kj (z) cos 



Kj iz) cos 



— Kj (z) cos 
9i(z 



4 

9i(z) 



9i (g) 
4 

9i (z) 
4 

COS^Q) 



sin On 



Here the Green's functions are given as, 
9 i}z + e) = 
9r(z) = V 



9r {z) + igi {z) , 
00 da; cosci;2; 



9i [z) 



-V 



-oo 27ri 1 - e-^ 
00 du sinu;z 



0J] 



2ml 



-0UJ 



9 (A) M 



( W ) _ g W {u) 



cos ^o(3.13a) 
cos^ (3.13b) 
(3.13c) 
(3.13d) 

(3.14a) 
(3.14b) 

(3.14c) 
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where g^ R \oj) = g(uo + iO + ) and g^ A \uj) = g(uj - iO + ). 



Fig. 5. Integration contour used when deriving eqs.(3.14) from eq.(3.13). 



In the next chapter, we will solve these equations and calculate the tunneling rate at 
absolute zero. However, since the integral of kernel Jq drK (r) in eq.(2.42) is logarithmically 
divergent at absolute zero, we will use finite temperature form for K{iz) in the calculation. 

3.2 Decay Rate Formula 

As far as the renormalized parameters are real, we can make use of the Korshnov's formula 
for the tunneling rate of the phase in the large dissipation limit©, and this condition is 
achieved in the neighborhood of the critical current. We find for the tunneling rate, T, 

T = /qmexp (S B ) , (3.15a) 
S B = -2a ren logsin# , (3.15b) 




Ax = (3.15d) 

V 1 - 167r2m r 2 en^ren/«r 2 en 

Here a ren , m ren , fl ren and 9 are determined from eq.(3.13) 

On the other hand, when the renormalized parameters acquire imaginary part, the free 
energy obtains imaginary part correspondingly. The appearance of the imaginary part re- 
flects the fact that the potential well is no longer metastable, but rather the phase can decay 
directly without the barrier. Hence we can regard this imaginary part of the free energy as 
the direct decay rate of the phase. The calculation of the imaginary part is performed by 
saddle point approximation. In this case the partition function will be estimated by changing 
the contour as in Fig. 6, 

Z = e-? F 
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Fig. 6. The path used in the saddle point approximation of eq.(3.17)(dot-dashed line). The saddle point 
is 6 - 



d9 



~ e 



2ir(3/m 
0w(e o ) 



-PW(6) 

iC/2 
-iC/2 



dy 



--e 2 



w"(e )i 



(3.16) 



'2irP/m 

where ( is the argument of W(8q). From this we obtain for the imaginary part of free energy 



as. 



ImF = ImW^o) + ~. 

P 



(3.17) 



Here ImW is given by following formulae ,0) 



ImW (e ) 

lmW l 

lmW 2 

ImW 3 
lmW 4 



Im [Wi + W 2 + W 3 + W 4 ] 



— Im 



log- 



m r 



m 



+ log 1 



7: 



K'H _|_ ^rcn 



Im 



— log 

2n —u< 



m 1 m 7 - 7 rcn - Q 
Im h — Im 



2 

ren 



2/T m T 

— — Im6 , 

2e 

-Imm reil fiL- 



27T m r 



LU + UJ C - UJ + U)_ 

— log — log 

Z7T — U) + Z7T 



— log — 

2tt & - 

- ^- 
— C<J_ 



Here w c is the cut-off frequency, 7 ren , o; + and u;_ are given by, 



Tren 

A, 



27rm r 

It 



\ 



[-1±A 2 ] 



^ ''ren 
/ ren 



(3.19a) 
(3.19b) 

(3.19c) 
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Decay rate is obtained by the relation T = 2lmF from eq.( 3.17|) . 



In the next section, we will apply these formula to the effective action, eq.(2.35), and 
investigate the decay rate of the phase in SNS junction. 

§4. Results and Discussion 

In this section, we will explain the results obtained from the effective action, eq.(2.35), by 
the approximation presented so far. For the numerical calculation, we have taken L = 0.4/im, 
W = 80/im and C = 200fF. 

4-1 Phase Diagram for j = 

We discuss the case j = in this subsection. A schematic representation of the phase 
diagram in Ej -a plane is shown in Fig. 7. The white region is for the localized phase(f2 ren 7^ 
0), where the phase variable is trapped in one of the local minima of the cosine potential. 
The hatched region, on the other hand, expresses the delocalized phase(f2 ren = 0). In the 
shaded region, our effective action fails, since we have an unphysical result, a ren < 0. The 
reason for the failure of our action in this region will be discussed in subsection 4.3. 



Fig. 7. The schematic representation of the phase diagram for j = plotted in the plane of Ejq and a. 
The hatched region is for the delocalized phase and white region for the localized phase. In the shaded 
region, our approximation fails. 



Two trivial limiting cases exist in eq.(3.13), i.e., Ejq = and a = 0. In both cases, we 
have K{iz) = 0, which leads to the unrenormalized result, m rGn = m, a ren = a and f2 ren = 0. 
This means the phase is delocalized. The delocalized phase extends to the finite area for 
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finite values of a and Ej . The boundary between the localized and delocalized region lies 
around a ~ 1, when Ej is small and comes closer to a = as Ej Q increases. This is due to 
the suppression of quantum fluctuation by the potential barrier. 

For larger a, the phase is localized in one of the minimum of the cosine well. In this region, 
the effect of cosine term is to increase m ren and to decrease a mn . In particular, on increasing 
Ejo, one observes that a ren decreases and finally becomes and eventually negative. We can 
not treat the region with larger E J0 .(See section 4.3.) 

We will show the a-dependence of renormalized parameters along lines A and B in the 
following. Those along line A are shown in Fig. 8. In this fig Ure, 77Z ren , Q^ren? ^ren and 
m ren ^ren are scaled with unrenormalized value, i.e., m, a, ojj and muj = Ej, respectively. 
Here Josephson frequency, uj, is defined as ujj = y/2E c Ej and Ej = Jq &tK{t). 



Fig. 8. Renormalized parameters m lcn /m, a rcn /a, fl rcn /u>j and m ren fl^. cn / Ej are plotted as a function of 
a along line A. 



The behavior observed in Fig. 8 can be understood as follows. In large a limit, K(iz) 
becomes short ranged, as was explained in chapter 2, and the retardation effect will not 
appear in the renormalization of m ren , a ren and f2 ren - Thus, in this limit, m ren , a ren and Q ren 
approach to m, a and u>j, respectively. On the other hand, in small a limit, the quantum 
fluctuation of the phase becomes large and it smears out the effect of cosine term, and we 
obtain m ren — > m, a ren — > a and f2 ren — ► 0. This leads to the localization-delocalization 
transition as a function of a, which was also found in a system of a particle moving in a 
periodic potential with Ohmic dissipationBil' For general values of a, the effect of cosine 
term is to increase m ren and decrease a ren . This tendency can be explained as follows. By 
considering the small fluctuation of the phase around 9 = 0, we can expand cosine term in 
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the action, eq.(2.30), as, 



Sj = - drdr'K(r - t') cos ^ 

~ - ^ drdT'K(r -r') 
Jo 



0(t) + 9{t')' 



'0(t) + 9(t')' 



— J29 n [K{iu; n = 0) + K(iu n )]9_ n , 

n 



(4.1) 



where we have neglected the constant term. By the analytic continuation, iu n — > u + i0 + , 
we find the correction to m and a by the kernel are given by the low energy behavior of 
Re K {R \uo) - KW(0) < 0, and ImK( R \u) > 0, respectively. The signs of these terms 
account for the enhancement of m ren and reduction of a ren - 

Next we will show in Fig. 9 the a dependence of the parameters on line B. For large a, 
the behavior observed on this line is essentially the same as that observed on line A. In this 
case, however, as a is decreased, we observe that a rcn becomes smaller and finally vanishes. 
For smaller a, we can not use our effective action. (See section 4.3.) 



Fig. 9. Renormalized parameters m len /m, a lcn /a 7 Q Ten /u)j and m ren fl^ en /Ej are plotted as a function of 
a along line B. 



4-2 Bias Current Dependence of Renormalized Parameters and Decay Rate 

In this section, we will introduce a finite bias current j and study the j-dependence of 
variational parameters and the decay rate of the phase from the metastable state. 



In Fig. 10, j dependences of renormalized mass, m ren , dissipation, a ren , frequency, Q r 
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Fig. 10. Rcnormalized parameters (a) m rcn /m, (b) a rcn /a, (c) fl len /u)j(l — (j /j c ) 2 ) 1//4 and (d) 
§0/ sin~ (j/j c ) plotted as a function of bias current j, for the points (a)~(d) of Fig. 7. 

and the position of the minimum, 6 , are shown. The parameters for these figures are 
Ej = 1.0K for all graphs, and a = 5 (a), a = 10 (b), a = 20 (c) and a = 120 (d). The 
parameters m ren , a ren , f2 ren and 60 are compared with the unrenormalized values, to, a, 
ujj(1 — (j/jc) 2 ) 1//4 and sin _1 (j/j c ), respectively. Here j c is the classical critical current given 
by j c = 2eEj. We notice that increasing the current, the renormalization of m ren and a rcn 
becomes smaller. This is because the factor cos^o i n eq.(3.13) decreases as j is increased 
and also the quantum fluctuation reduces the effect of cosine term. Therefore, m ren and 
a ren approach to their classical value. On the other hand, quantum fluctuation reduces f2 ren 
compared to its classical value, but enhances 8q. When the current is increased above a 
critical value, jf^, no real physical solution to the self-consistent equation exists and the 
parameters acquire imaginary part. This means that the localized state in the well is no 
longer metastable but becomes unstable due to zero point fluctuation and decays directly 
from the well. In the following, we will call this process as direct decay and its onset jf 1 as 
effective critical current. 

Now let us discuss the bias current dependences of the effective critical current, j^ s , and 
the decay rate, T, with a and E J0 fixed. 
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To give a reference frame, we will compare the result with the approach where no re- 
tardation effect is taken into account. One is the local and semi-classical approximation, 
in which the kernel K{r) is treated as a 5-function and the fluctuation of the phase is 
neglected. (Denoted SC in the following.) Second is the local approximation in which the 
quantum fluctuation are included by self-consistent harmonic approximation as was done 
by Kleinert, but the retardation effect of the kernel was neglected. (Denoted 'without retar- 
dation' in the following.) Here the fluctuation around the self-consistent solution, which 
was taken into account in the Kleinert's original approach, is not considered. (See section 
4.3.) The last is our approximation, in which both the retardation effect of the kernel and 
the quantum fluctuation are taken into account. We fix the classical Josephson energy, 
Ej = Jq K{r)dr, for the above comparison. 

It is to be noted that these three different theories treat the renormalization of the param- 
eters due to the kernel, K(r), and quantum fluctuation very differently. These renormaliza- 
tion are not considered in the first approach(SC), while in the second treatment ('without 
retardation'), there is no renormalization to mass and dissipation but the renormalization 
of frequency, since the kernel has no time dependence, and only Q rcn and # are treated as 
variational parameters. 

In Fig. 11 and Fig. 12, j dependence of the decay rate T is shown for the points (a) and 
(d), respectively. (See also Fig. 7) Here, T is plotted in unit of uj. Since the parameter a is 
not large enough to use Korshnov's formula, we have used formula for cubic plus quadratic 
potential for the tunneling rate in Fig. 11. HH) The solid line is our result, dotted and the 
dot-dashed lines are by SC and 'without retardation', respectively. 



Fig. 11. Decay rate for point (a) in Fig. 7 plotted as a function of j/j c - 
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Fig. 12. Decay rate for point (d) in Fig. 7 plotted as a function of j. 



The overall behavior of the decay rate of our result is similar to that of 'without retarda- 
tion'. With the increase in j toward jf 1 , decay rate increases. For current larger than jf 1 , 
there is no contribution from the bounce solution. Instead the contribution from another 
process, i. e. the direct decay, sets in. The direct decay rate is far smaller than the tunneling 
rate. This may partly be due to the overestimation of the tunneling rate in the localized 
regime as will be discussed later. 

Next we discuss the difference between 'without retardation' and our results in jf* and F" in 
detail. The j-dependence of decay rate for the parameters corresponding to the parameters 
of (d) in Fig. 7 is shown in Fig. 12. Here we notice that jf 1 is larger in our case than in 
'without retardation'. This implies the smaller quantum fluctuation for non-local kernel 
for point (d). As a result T is smaller in our case than in 'without retardation'. On the 
other hand, at point (a), we observe the opposite behavior for jf 1 and r.(See Fig. 11.) The 
different behavior seen in the quantum fluctuation of points (a) and (d) can be considered 
as a result of competition between the effect of m rcn and a rcn . The renormalization effect 
makes m ren larger and a ren smaller. The former suppresses the quantum fluctuation and the 
latter enhances it. For large a, the former effect is dominant, while for small a the latter 
effect becomes stronger and jf 1 of our result is suppressed compared to that of 'without 
retardation'. 

The a dependence of j^/jc along line A is shown in Fig. 13. The white circles represent 
our result and the black ones for 'without retardation'. We observe that the suppression of 
the critical current is not so large, except in the neighborhood of localization-delocalization 
transition. The dependence of jf 1 on the parameter, a, is qualitatively the same as the 
classical one. The difference in jf 1 between our results and 'without retardation' is also 
small. 
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Fig. 13. The a dependence of the effective critical current along line A. White circles denote our result 
and black ones the result of 'without retardation'. 



4-3 Discussion 

We conclude that the effect of the retardation of the kernel, K(t), on the renormalization 
of jf 1 is small compared to the case of local kernel K(r) = Ej5(t). The most significant 
effect, however, appears on the renormalization of mass and dissipation. The retardation 
effect leads to the j dependent renormalization of mass and dissipation, which does not 
appear for the local kernel. 

Next we will discuss the limitation of our effective action. As shown in the phase diagram 
in Fig. 7, we observe the renormalized dissipation, a rcn , becomes negative in the shaded 
region. Since we have derived the effective action in the lowest order perturbation theory 
with respect to the matrix element of transfer integral of the electron at the SN boundary, 
the parameters a is given by the inverse resistance, 1/Rn, of N region without coupling to 
the superconductor. We expect the terms which compensates the negative contribution from 
the cosine term appears, when the higher order terms in perturbation series are considered. 
Thus, we believe that the negative dissipation is an artifact of our derivation of effective 
action, in which only the lowest order terms are considered. 

Finally, we comment on the effect of fluctuation around the self-consistent solution. In 
our approach, tunneling rate, T, is evaluated with the renormalized potential obtained in 
the self-consistent harmonic approximation. When the effect of fluctuation around the self- 
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consistent solution is considered, the tunneling rate is expected to be reduced. It is difficult, 
however, to treat this effect in our effective action, eq.(2.35). Instead we here consider the 
simple case of cubic potential 'without retardation', i.e., 

V(x) = -mwlx 2 {\ - — ). (4.2) 
2 xo 

When j ~ j c , this gives a good approximation to the cosine potential. The renormalized 
potential, V(x), is given by, 

1 1 uo 2 

V(x) = -mfl r 2 m i 2 - -m^x 3 , (4.3) 
Z 2 xq 

where the effective attempt frequency, f2 ren , is smaller than uj$. The imaginary part of free 
energy is calculated from the partition function, 

Z = I Vxe~ s 

Vxe~ s e- {s - s \ (4.4) 



where S is given by replacing V(x) in S with V(x). In the calculation of tunneling rate, we 
have evaluated S by the bounce solution. Taking the fluctuation into account, we observe 
that additional factor, e~ ( ~ s ~ s ' }B , arises in the tunneling rate formula. Here (•■ -}b denotes 
the expectation value with respect to the fluctuation around the bounce path without the 
zero energy fluctuation. Since the expectation value (y 2 ) b is positive, the expectation value 
of the factor, 

S - § = l -m{ul - n 2 J f dr(x B + y) 2 , (4.5) 



2 -\— u reny i 

JO 

is positive. Thus the tunneling rate is suppressed, when we consider the fluctuation around 
the self consistent solution. This situation is the same for 4 potential©. We also expect 
that this situation does not change when we consider the effective action for SNS junction 
and the tunneling rate is suppressed to be closer to the direct decay rate. 

§5. Summary 

In this paper, we have investigated the effect of low energy excitation due to conduction 
electron on the dynamics of the phase in SNS junction. We derived the effective action for 
the phase in SNS junction from a microscopic Hamiltonian. The action is different from that 
of Josephson junction in two respects. First, there is a term describing Ohmic dissipation 
expressed by the normal resistance of the junction, which is the resistance of the junction for 
large bias current limit. Second, the kernel, K(t), which describes the Josephson coupling, 
has a long time tail. Both reflect the existence of low energy excitation in N region. The 
adiabatic approximation cannot be justified a priori to this action, and we focused on the 
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effect of retardation of K{r) on the dynamics of the phase. To clarify the effect, we have 
calculated the tunneling rate of the phase out of a minimum of potential with and without 
the retardation effect of kernel. 

The case where the N region is formed by the diffusive two dimensional electron gas was 
examined from the experimental interest. The time dependence of K(t) is characterized by 
the diffusion constant, D, which is proportional to the strength of the Ohmic dissipation a. 
Writing K{t) as K{r) = Ej k(r), where Ej Q is a multiplication factor and k(r) is a function 
of a, we have investigated dynamics of the phase by changing E J0 and a. 

The effect of the kernel was examined in the self-consistent harmonic approximation, and 
was expressed by renormalization of mass, m ren , dissipation, a mn and the attempt frequency, 
fl r en, of the phase variable. We have calculated the tunneling rate, T, of the phase by 
these renormalized parameters on the plane of a and Ejq. These results are compared 
with the previous approaches; one is the semi-classical approximation(SC) in which neither 
of the retardation effect and quantum fluctuation is considered, and the other is the local 
approximation (without retardation), in which only the quantum fluctuation is considered. 

In large a limit, the motion of the phase is classical and our result approach to the semi- 
classical calculation. On the other hand for a < 1, the quantum fluctuation smears out the 
effect of cosine term. For the intermediate region, we found bias dependent renormaliza- 
tion of dissipation and mass. The renormalization of dissipation makes the tunneling rate 
larger and that of mass makes it smaller compared to the unrenormalized case resulting in 
delicate competition of these two factors. For a large value of dissipation, the latter effect 
dominates and the retardation effect makes the tunneling rate smaller. For small value of 
the dissipation, vice versa. 

In conclusion, we have investigated the effect of retardation of the kernel K(t) on the 
dynamics of the phase. We find the effect appears most significantly on the renormalization 
of mass and dissipation. 

Current dependent renormalization effect can be detected in the experiments of AC Joseph- 
son effect by measuring Q-value of plasma resonance directly. It is possible to measure 
plasma frequency u> p and the Q-value u p RC at the same time by applying small power rf 
voltage to the junction as in Dahm et al&. Since C and R are proportional to m ren and 
1 1 a respectively, it is expected to observe the reduction of RC value as the increase of bias 
current in SNS junction. As for the renormalization of j% and tunneling rate, T, however, 
the dependence to the time scale of the K(r) is found to be small. 
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Appendix: Analytic Continuation of the Self Consistent Equation 

In this appendix, we will show the procedure to derive eqs.(3.13) from eq.(3.12) by ana- 
lytical continuation. 
The self consistent equation in imaginary time is given as, 
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where the function g(r) and K(r) are given by, 
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Our aim is analytic continuation of eq.(A-l) from imaginary frequency to real frequency. 
We first change the time variables as iz = r. With this, £ (Lv n ), becomes, 



£ (ico n ) = - drK (r) (l + e^A e -i^(o)+ s (r)} 
2 Jo v ' 



\dzK (iz) (1 + e WnZ ) e-i^(o)+ ff (i^)}_ 



(A-3) 



Since there is no singularities in the region surrounded by C, C_, Coo and Cp(Fig.5), we 
can change the integration contour as:C — > C_ + + Cp and obtain, 
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Here the contribution from is zero. The analytically continued form of the functions, 
K(iz) and g(iz), are given by, 



K (iz + e) 
K R (z) 
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Note that when all the renormalized parameters are real, the definition g R and gj corresponds 
to the real part and imaginary part, respectively. Otherwise, these functions have imaginary 
part. Since \voK^(uj) is an odd function of uj, we can show following relation easily, 



K R (z) = 

Kj(z) =- 
K(iz + /3-e) = 
K (iz - e) = 

The same relations as eq.(A-6) hold for g(iz). 
With eqs.(A-5) and (A-6), eq.(A-4) becomes, 
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By the analytical continuation, io; n — > uj + iO + , we obtain, 
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At low frequency, we can expand eq.( [A-8| ), and obtain the self-consistent equations, 
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Hence eq.(3.13) is derived from eq.(3.12). 
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